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Abstract 

We consider the possibility that photons of noncommutative QED can make 
bound states. Using the potential model, developed based on the constituent gluon 
picture of QCD glue-balls, arguments are presented in favor of existence of these 
bound states. The basic ingredient of potential model is that the self-interacting 
massless gauge particles may get mass by inclusion non-perturbative effects. 
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1 Introduction 


In Abelian gauge theories on ordinary space-time, there is no self-coupling between the 
gauge fields. The best known example is the quantum theory of interaction between 
electric charges and photons, the so-called QED. The situation is different in non-Abelian 
theories, and due to the commutator term in the held strength, these theories are involved 
by direct interaction between the quanta of gauge helds. It is now widely believed that the 
strong interaction is described by a non-Abelian gauge theory accompanied with proper 
matter helds as quarks, the so-called QCD. As gluons are the quanta of QCD gauge 
held, from the very beginning this possibility was considered that if gluons can make 
bound states free from valance quarks, the so-called glue-balls. Although the properties 
of glue-balls have been studied for a long time, their existence have not still been approved 
experimentally. 

Recently a great interest has been appeared to study held theories on spaces whose 
coordinates do not commute. These spaces, as well as the held theories dehned on them, 
are known under the names of noncommutative spaces and theories. In contrast to QED 
on ordinary space-time, as we briehy review in next section, noncommutative QED is 
involved by direct interactions between photons. Interestingly one hnds the situation 
very reminiscent to that of non-Abelian gauge theories, and then the question is whether 
there are some kinds of bound states in analogy with glue-balls of QCD, here might be 
called “photo-ball” s. It is this question that we consider in this work. Our approach to 
study photo-balls is based on one of methods that has been developed for glue-balls. As 
glue-balls are non-perturbative in nature, there is still no systematic way for calculation 
of their properties from the hrst principles of QCD. Instead, among the years many 
approaches are developed for extracting glue-ball’s properties, though each approach is 
based on expectations or estimate calculations. 

Among many others, one approach for studying properties of glue-balls has been the 
so-called constituent gluon model. In any study of bound state of gluons, one is encoun¬ 
tered with a situation in which gluons, though at hrst were introduced massless to the 
Lagrangian, are bound and do not disjoint to propagate to inhnity. Correspondingly, it 
is argued that quantum huctuations around a charged particle, that should be treated 
non-perturbatively in QCD, can make an accompanied cloud for it, causing a dynamically 
generated mass DEa Accordingly, it appears very useful to dehne constituent quarks 
and gluons, for which we assume a mass of order of bound states of the theory, while their 
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Lagrangian counterparts may be massless or almost massless. As extracting the masses 
of constituent particles from first principles is not yet done in a satisfactory way, the best 
evaluations come from estimations based on general considerations, phenomenology or 
lattice calculations. 

Once one accepts that a glue-ball is a bound state of constituent gluons, the question is 
about the effective theory that captures the interaction between them. One approach is to 
consider constituent gluons as massive quanta of an effective gauge theory. It needs some 
kinds of proof, but hopefully this effective gauge theory has the same qualitative features 
of the true (massless) theory, but in non-perturbative regime Pl E] . Since it is believed 
that the main contribution to the mass of a glue-ball is coming from the constituent 
masses of gluons, it is expected that constituent gluons move non-relativistically inside 
the glue-ball, and so perturbative calculations for hnding the effective potential should be 
done in non-relativistic regime. Having the effective potential at hand, by studying the 
Schrodinger-type equations, one can have estimations about mass or size of glue-balls. It 
is the heart of the potential model approach for studying the properties of the glue-balls 

iHiiiin]. 

There are two related issues when we are considering the effective gauge theory of 
constituent gluons as massive vector particles. First, it is known that the gauge symme¬ 
try is lost via the mass term, and the second, massive gauge theories are known to be 
perturbatively non-renormalizable. Here we give comments on these issues ^ E] . The 
non-renormalizability of massive gauge theories is under this assumption that the mass 
in the theory appears as a hxed parameter, surviving at large momentum. In fact the 
insufficient decrease of propagator of a massive vector particle at large momentum, due 
to simple power counting, suggests that the theory can not be renormalizable. But the 
situation might be different in a theory with constituent mass. At very large momentum, 
where coupling constant is small due to asymptotic freedom, the perturbation is valid and 
gluons appear as massless particles. So the mass of constituent gluon, which is generated 
dynamically, depends on momentum and vanishes at large momentum. In a theory for 
gluons, it is argued that if one can keep the dependence of constituent mass on momen¬ 
tum, which of course is possible only by including the non-perturbative effects, the theory 
may appear to be non-perturbatively renormalizable. 

Although the argument above is for a model involved by dynamically generated mass, 
due to lack of a systematic treatment of non-perturbative effects, much can be learned 
via a kinematical description of gluon mass |2]; it is to assume mass as a £x parameter. 
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though the problem still remains with local gauge symmetry. To overcome this problem, 
there is a prescription that we review briefly in below. Let us consider firstly the case with 
Abelian symmetry, and QCD case appears just as a generalization IIII2I. The starting 
point is the theory given by Lagrangian density: 



in which g is the coupling constant, and is a scalar held. The Lagrangian is invariant 
under local symmetry 

A'^ = A^- (p' = (p - A, (2) 

with A as symmetry transformation parameter. Now we see although the gauge held 
has got mass, the local symmetry is kept. Of course we mention that giving mass is 
done with the price of introducing an extra scalar held. We have another example of 
this observation in spontaneous symmetry breaking mechanism, in which we are left 
with Goldstone bosons. In fact, this extra scalar held, just like their Goldstone boson 
counterparts, does not appear in the A-matrix, i.e. as external legs of diagrams. One 
may insert the scalar held into the Lagrangian via the obvious solution: 

(p = g^d-A, (3) 

getting a non-local but still gauge invariant theory involving only A^. It is the general¬ 
ization of this mechanism that has been used for the case of constituent gluon description 
of QGD glue-balls ISHUE], and here we use for photo-balls of noncommutative QED. The 
starting point for the QGD case is the Lagrangian density 1110: 

L = (4) 

in which 


V (cp) = exp ^ TV, 

a 


( 5 ) 


and F/j,iy = dfj^A^ — d^Afj, — ig[A^, , with [T“, T^] = The action is invariant under 

- 4 ;. = uA^u-^ + tua^u-\ 


V{^') = UV{^), 


( 6 ) 
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in which U = U{x) is unitary matrix defining the transformation. Again one can insert 
the scalar fields into the Lagrangian via power series solution in g 1110 : 

= 3^9-A-3^l---U- (7) 

As mentioned above the extra scalars do not appear as external legs of diagrams, but the 
situation is even simpler as far as one considers just the tree diagrams, in which one can 
ignore the scalars. So for tree diagrams, and in proper gauge, the Lagrangian density in 
use is practically mu- 

L = ( 8 ) 

simply as a gauge theory for massive gluons. 

In the non-relativistic limit the potential can be read off from the total invariant 
amplitude M.fi via the Fourier transform jSlIllEl: 


l/(r) 


d^q 

—r — , = iM. fi 

87r3 A^EiE2E^Ei ^ 


(9) 


in which q is the momentum transferred between the in-coming particles. The total 
invariant amplitude gets contribution at tree level from the s-, t-, u- channels, and the 
so-called seagull (s.g.) diagram, coming from 4-gluon vertex of QCD 0111 El- In non- 
relativistic limit it can be shown that the s-channel’s contribution is ignorable, to get the 
final expression 000: 


• O J^nr'p JjhHp 

'‘Mf, = (4m= + 3q"-2SV + 2(S.q)= + 6»S.(qxp.: 


q2 _|_ ^2 

- 2)j. 


( 10 ) 


By above the gluon-gluon potential in one-gluon-exchange (oge) approximation appears 
in the form 000 : 


h^ge(r) — 


_ g2 j^ace jbde 

dvr 


1 lo2 
- -h -S^ 
14 3 



^(L.S)V-i™=((S.V) 


2 £abe £cde 


g^r^f 


IT 


dvr 




1 

3 




e 


r 


( 11 ) 


in which Vr = r ^dr- 
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The organization of the rest of this work is as follows. In Sec. 2 we review some basic 
features of canonical noncommutative spaces, and also the held theories dehned on them, 
specially noncommutative QED. We also remark on some aspects of noncommutative QED 
that make this theory in some extents similar to QCD. Sec. 3 is devoted to extracting 
the effective potential between massive photons, by studying the non-relativistic behavior 
of their scattering. Sec. 4 mainly contains the dynamics of photons under the effective 
potential obtained in Sec. 3. The existence proof of bound states is also presented in Sec. 
4. Sec. 5 is for conclusion and discussions. 


2 Noncommutative Space-Time And QED 

In last years a great attention is appeared in formulation and studying held theories 
on noncommutative spaces. One of the original motivations has been to get hnite held 
theories via the intrinsic regularizations which are encoded in some of noncommutative 
spaces jH] . The other motivation comes back to the natural appearance of noncommutative 
spaces in some areas of physics, and the recent one in string theory. It has been understood 
that string theory is involved by some kinds of noncommutativities; two examples are, 1) 
the coordinates of bound-states of N D-branes are presented by iV x iV Hermitian matrices 
[Zj, and 2) the longitudinal directions of D-branes in the presence of B-held background 
appear to be noncommutative, as are seen by the ends of open strings |H1 El CD]- In 
the latter case for constant background one simply gets the canonical noncommutative 
space-time, introduced by commutation relations for coordinates as: 

( 12 ) 

in which the parameter a; has the dimension of energy, and signihes the scale where 
noncommutative effects become relevant. The is a real antisymmetric matrix with 
elements of order one. Now since the coordinates do not commute, any definition of 
functions or fields should be performed under a prescription for ordering of coordinates, 
and one choice can be the symmetric one, the so-called Weyl ordering. To any function 
/(x) on ordinary space, one can assign an operator 0/ by 

6/(4 /«/>)<*■“ ( 13 ) 
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( 14 ) 


in which f{k) is the Fourier transform of f{x) defined by 

f(k) = j d"x f{x) 

Due to presence of the phase in definition of O/, we recover the Weyl prescription 

for the coordinates. In a reverse way we also can assign to any symmetrized operator a 
function or field living on the noncommutative plane. Also, we can assign to product of 
any two operators Of and Og another operator as 

Of ■ Og =: Of^g (15) 


in which / and g are multiplied under the so-called ^-product defined by 

iQfJ'k' 

= exp {^d^^dy^)f{x)g{y) \y=^ (16) 

By these all one learns how to define physical theories on noncommutative space-time, and 
eventually it appears that the noncommutative field theories are defined by actions that 
are essentially the same as in ordinary space-time, with the exception that the products 
between helds are replaced by ^-product; see m as review. Though ^-product itself is 
not commutative (he., fg g ^ f) the following identities make some of calculations 
easier: 


f*gd-x= / g^fd-x= / fgd-x 


f-kg-khd"^x= / f{g-kh) d'^x= {f-k g)h d”x 


f -k g -kh d"'x = / h-k f k g d”x = g k hk f d"'x 


By the first two ones we see that, in integrands always one of the stars can be removed. 
Besides it can be seen that the T^-product is associative, he., fkgkh = {fkg)kh = fk{gkh), 
and so it is not important which two ones should be multiplied firstly. 

The pure gauge held sector of noncommutative QED is dehned by the action 

Sg.ug.-a.M = -ly d'l *F'“' 

= -i y d-i (17) 

in which the held strength Fgu is 


= dyA^{x) - d^Ay{x) - ie[Ay{x), A^{x)]^ 


(18) 


7 



by definition [f,g]* = f ^ g — g ^ f ■ We mention [x^,x^]i, = iO^'". The action above is 
invariant nnder local gange symmetry transformations 

A'^{x) = Ui<A^{x)i.U-^+ -Ui.d^U-^ (19) 

in which U = U{x) is the *-phase, defined by a fnnction A(a;) via the ^-exponential: 

U{x) = exp^(iA) = l-fiA — - AtAtA-I--- - , (20) 

= 1 ( 21 ) 

in which U~^ = exp^(—iA). Under above transformation, the field strength transforms as 

= ( 22 ) 

We mention that the transformations of gange field as well as the field strength look like to 
those of non-Abelian gange theories. Besides we see that the action contains terms which 
are responsible for interaction between the gange particles, again as the sitnation we have 
in non-Abelian gange theories. We see how the noncommntativity of coordinates indnces 
properties on fields and their transformations, as if they were belong to a non-Abelian 
theory; the snbject that how the characters of coordinates and fields may be related to 
each other is discnssed in ini. These observations make it reasonable to stndy whether 
and how the photons can make bound states in such a theory. 

There is another observation that promotes the formal similarities of noncommutative 
and non-Abelian theories to their behaviors, that is the negative sign of /3-function, which 
manifests that these theories are asymptotically free [ISIII3. By this it is more reasonable 
to see if the techniques developed for QCD purposes can also be used for noncommutative 
QED. 

The phenomenological implications of possible noncommutative coordinates have been 
the subject of a very large number of research publications in last years. Among many 
others, here we can give just a brief list of works, and specially those concerning the 
phenomenological implications of noncommutative QED. The effect of noncommntativity 
of space-time is studied for possible modifications that may appear in high energy scat¬ 
tering amplitudes of particles ra. in energy levels of light atoms HHnni, and anomalous 
magnetic moment of electron QSI. The ultra-high energy scattering of massless photons 
of noncommutative QED is considered in US and the tiny change in the total amplitude 
is obtained as a function of the total energy. Some other interesting features of noncom¬ 
mutative ED and QED are discussed in [20]. The issue of formation of new bound states 
by space-time noncommntativity has been considered in EH- 


3 Massive Noncommutative QED And Effective 
Potential Between Photons 


3.1 Massive Photon-Photon Scattering Amplitude 

Now here, following the procedure developed for QCD case, we give mass to photons 
of noncommutative QED. As described this is done by introducing an extra scalar field, 
getting the Lagrangian density: 

-^ = -^gauge-field H-H- d^V{(p)'kV (23) 

in which (■■■)* = (■■■)*(■■■), and V{(p) is the T*r-phase dehned by the scalar held ip; see 
(EDI). The action dehned by above Lagrangian is invariant under transformations 

V{ip') = Ui.V{ip) (24) 


in which U is the same in (HI. Now we just list the Feynman rules mniiiHi We choose 
the gauge in which the propagator takes the form: 






+ 




p2 _ ^2 (jp _ 

In the non-relativistic limit we have for momentum and polarization vectors iis] 


(25) 


A + I^.P) 


P) 


m 2m? 

in which e is a 3-vector satisfying e* ■ e = 1. From Lorentz condition laHHIHI, we have 


(26) 

(27) 


p-e = p^e^ = 0 (28) 

In this work we assume for the signature of metric = (+1, —1, —1, —1). As in this work 
we restrict ourselves to tree diagrams, after removing one *, and by Lorentz condition, 
we practically are using the Lagrangian [2] 

L = (29) 
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Figure 1: 3- and 4- photon vertices of noncommutative QED. 


with dehned in m- There are 3- and 4- photon vertices given in Fig. 1. As in this 
work we consider noncommutativity just in spatial directions, that is assuming 6^^ = 0 
for i = 1,2, 3, we have for the vertex-functions [IHIIH!: 

Ici [>^ Ico r 

= -2esin(^^) {h - + {h - + {h - (30) 


and 


^ kiMMM ~ 


. .ki x k 2 x . /ka X k 4 . , 
sm(---) sm(---) [g 


+ sin( 


2 

ka X ki 


) sin( 


2 

k2 X k4 




)(9 


sm 


2 ' ' 2 
ki X k4. . .k2 X ka. . 

-3- 


M1/44 ^/42 M3 _ ^ A* 1 M2 ^ M3 M4 j 


MlM2^M3M4_^MlM3^M2M4 


(31) 


in which a x b = 9'^^aibj, and the momenta and indices are given in Fig. 1. Also in each 
vertex energy-momentum conservation should be understood. 

Although there are four diagrams at tree level, those coming from s-, t-, u- channels, 
and seagull diagram (Fig. 2), when extracting the potential, by the properly symmetrized 
wave function for identical particle systems, the “exchange” or “symmetry” diagrams are 
automatically taken care of, causing that only one of t- and u- channels’ contributions 
should be added to others’ contributions jl]. For s-channel we have 


iM)i = 4e^sin( ^^ ^ sin( P^^ ^ 


~'^gpLv 




2 

+ 


ik^kjy 


m 


pJ/M3M4 


2 ^ (P -fc,P3,P4 3/.3 4^4 


(32) 
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s-channel 


t-channel 


seagull 


Figure 2: s-channel, t-channel, and seagull diagrams. 


In s-channel we have k = pi + p 2 = Ps + P 4 , and by using it we get 

• i /P3i<k\ 

= -Aie ^— % _ ^ ' ^ 2(^1 - P2)p + 2e2pP2 ■ ei - 2eipPi ■ 63 ) 




(^63 ■ 64(^4 - V?iY - 2 e*Yp 4 ■ + 2 e*Yp 3 ■ elj 


(33) 


From the Lorentz condition the last term in propagator is omitted. By replacing the 
non-relativistic limit of e’s, we see that, even without considering coefficient involving 
sin(- ■ ■), the leading order contribution of s-channel is order of |pp m^, that we ignore 
comparing the zeroth orders. This observation is exactly the same that happens in the 
QCD case [31lll|3j. 

Now we come to t-channel. 


zA^}^=- 4 e^sin( ^^ ^ ^ )sin( ^^ " Ppxipi + P3)p + gxpipi - 2^3)^ + Pppip^ - 2pi)x 


P2 X q 


*A 

^1^3 ■ 


2 ' ' 2 


q2 _ ^2 


6364 


gua{P2 + PYs + ga6{P2 “ 2 ^ 4 ) 1 . + gSuiPi - 2p2)a (34) 


in which q = Pz — Pi = P 2 — Pa- It appears to be useful if we dehne in IS] 
(3| Jp|l) = 63 ■ ei(pi +P3)p - 2elpPz ■ ei - 2eipPi ■ 
and a similar one for (4|7*^12), and since (3| = (4|7^12)^^ = 0, thus, 

dn(El^)sin(EfS) 




q2 _ ^2 


-{3|JJ1){4|J0|2) 


(35) 


(36) 
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We continue in the center-of-mass frame, for which 


Pi = -P2 = Pi, 

P3 = -P4 = Pi + q = P/, 

q = P3-pi = P/-Pi- (37) 

In the final expression we keep just the zeroth order of |p|, although we should be careful 
about q dependence, on which we integrate over to find the large-distance behavior of 
the effective potential. So in the following steps we still should keep orders of |pp in 
{3|J,|l){4|J<'|2). 

The next steps of calculations are essentially those done for glue-balls that 

we present in below. Photon’s spin is one, and so it is useful to define the operators 
S = {S\S\S^) as 


= 

/ 0 0 0 \ 

0 0 -i 

, = 

/ 0 0 z \ 

0 0 0 

, 53 = 

/ 0 -i 0 \ 

i 0 0 


0 i 0 y 


—i 0 0 y 


\ 0 0 0 J 


for which we have [S\ 

Now we calculate the components of (3| Using (0!)|l and (ITTjl we get 

(3|J°|1) = -(2m + ^i^^)(e;-ei) 

+ —[-2(pi ■ ei)(p 3 ■ e^) + (p 3 ■ ei)(p 3 ■ e*) + (pi ■ ei)(pi ■ e^)] (39) 

TTT/ 

For any two vectors A and B we mention 

/l.B, = [A.B 1 -(S-B)(S-A)]_^ (40) 

(in the following we drop the identity operator 1). Using above we can write 

(3|J°|1) = -e^2m+ 

+—egfq^ - (S ■ q)^ + (S ■ pi)(S ■ Pa) - (S ■ P3)(S ■ pi)]ei (41) 

in which ci and 63 are column-matrix representation of 3-vectors ei and 63 . Again using 
the identity 

(S ■ pi)(S ■ Pa) - (S ■ P 3 )(S ■ pi) = p\pi[S\ 5^] = i{pi x pa) ■ S (42) 
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we finally reach to 


(3| J°|l) = -e^ {-2m- ~ ^ Ps) ' ^i] (^3) 

in which Si as hrst photon’s spin-operator. Similarly, 

(4| J°|2) = -el (-2m- ^ [q 2 _ (S 2 ■ q)^ -f f(p 2 x P 4 ) ■ S 2 ] 62 (44) 

- ZiTDj 7Y}j 

Now using p ■ e = 0, we can write the last two terms of (3| J^|l) in terms of q, and so we 
get 

(3|J^|1) = -e|(pi-h P 3 )^ei-h 2efq-Cl - 2eiq -63 

= -elipi + P 3 )''ei + 2el^[5^^q^ - 5^^q^]e{ (45) 

Let us dehne the matrices Bi, B 2 , B^ by their elements, 

B^ = 6^^q^ - P^q^ (46) 

and so for example we see for Bi that 

/ 0 q^ q^ \ 

fii = -g2 0 0 =-i[(SiXq)]^ 

y —q^ 0 0 y 

yielding 5'^^q^ — 5^^q^ = —i[(Si x q)]^^- By using the above we get 

(3|J|1) = e^[-(pi -h Ps) - 2f(Si X q)]ei 
(4|J|2) = ei[-(p2 + P 4 ) + 2f(S 2 X q)]e2 

By these all we have 

(3|J^|l)(4|Jp|2) = e^el[4m^ + p^ + p^ + p^ + p^-4q2 

-(Pi + Ps) ■ (P 2 + P 4 ) + 2(Si ■ q)^ -h 2(82 ■ q)^ 

-4(Si X q) ■ (S 2 X q) -h 2 i(p 3 X pi) ■ Si 2 i(p 4 X P 2 ) ■ S 2 
+2i{j)i + P 3 ) ■ (S 2 X q) - 2 i(p 2 -h P 4 ) ■ (Si X q)]eie 2 (50) 

in which we should care the order in multiplying the objects in above. By conservation of 
energy we have Pi + P 2 = pI + pI, and so in the center-of-mass frame Pi = P 2 = P 3 = pI- 


(47) 


(48) 

(49) 
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Also since (pi — P 3 )(pi + Ps) = 0 and so (pi + Pi — q) ■ (pi — Pi + q) = 0 we obtain 
Pi ■ q = ^q^ = -Pf ■ q- By using these all we have (pi + ps) ■ (p 2 + P 4 ) = -4pi + q^, 
and so 

P? + P 2 + Ps + P 4 - 4q^ - (Pi + Ps) ■ (P 2 + P 4 ) = 8p- - 5q2 (51) 


that we ignore pf in comparison with For the terms involving spin-operators defining 
the total spin S = Si -|- S 2 , we get Si ■ S 2 = |S^ — 2, and thus 


(Si ■ q)2 + (S 2 ■ q)' - 2(Si X q) ■ (S 2 X q) = (4 - + (S ■ q)' (52) 


By two other relations 


^(Ps X pi) ■ Si -h ^(P4 X P 2 ) ■ S 2 = i(q X Pi) ■ S (53) 

*(Pi + Ps) ■ (S 2 X q) - i(p 2 P 4 ) ■ (Si X q) = 2iS ■ (q X pi) (54) 

and using = 0 — q^ and Pi = — P 2 = P, we obtain 

iM\i = 4ie^ — J [4m^ -|- 3q^ — 2S^q^ -|- 2(S ■ q)^ -|- 6iS ■ (q x p)l -|- O(p^) (55) 

We mention that the kinematical dependence of t-channel amplitude, given by terms in 
[••■], no surprisingly is exactly that for gluons, presented in relation m- In fact the 
only difference between the case of gluons and photons in noncommutative QED is in the 
pre-factor, originated from difference between structure constants of group that appear 
in vertex functions m 

Now we come to seagull diagram, with the contribution 




r 

. /Pi X ; 

P 2 ^ 

. /Ps X 1 

P 4 n 


sm - 

^ 2 

) sin( 2 

) 


. /Ps X 

Pin 

. /P2 X 

P 4 

+ 

sm - 

^ 2 

) 

sm - 

^ 2 



. /Pi X 

P 4 n 

. /P2 X 

Ps 

+ 

sm - 

^ 2 

) 

sm - 

^ 2 



)(9'‘V”-9'‘'’9‘’t 

)( 9 '‘V"- 9 '‘V‘')] 


leading to 




= —8ie^[sin( 


Pi X P2 


) sin( 


Ps X P4 


-|- sin( 


2 ' ' 2 
Pi X Psx . /P2 X P4 


)ei 


. /J 

-j sm(- 


(56) 


) X (ei • €263 • €4 — Cl • €462 ■ 63)] ( 57 ) 
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If we keep only zeroth orders of momentum in the bracket, we have 



We can write as ei ■ 62 63 ■ 64 = Expressing and having 

from the spin-form {S^S’‘)ij = 5^^5ij — 5\5^, we get efe^ = e\{5^^ — S'lS{)ei. Using this we 
obtain 



In above, we use the relation S x S = iS. Finally we have 



(60) 


Now we mention that the contribution of seagull channel, for small noncommutativity 
parameter is something proportional to ( 6 'p)^p^ which is order of p^ that we ignore. 
This observation is different from that for QCD glue-balls, for them the contribution of 
seagull diagram is in zeroth order of momentum, and thus should be kept. The seagull’s 
contribution appears to be in form of 5(r) in the potential (ITT|l . 

3.2 Effective Potential Between Photons 

Before we proceed, we define the vector 6 based on tensor 6^^ by 



1 


(61) 


By this vector we can write the x-pro duct as 


a X b = = a*- eiij6’‘lP = - 0 ■ (a x b) = - b ■ (0 x a) 


(62) 


By this we have for the t-channel contribution: 



(63) 


in which A = | p x 0 and 

T(q) = 4m^ -|- 3q^ — 2S^q^ -|- 2(S ■ q)^ -|- 6*8 ■ (q x p) 


(64) 
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By the total amplitude the potential can be deduced using Q 


V2,{r) = 


d^q 


;p*qr 


le 


A\/E1E2ESE, 


AMfi 


‘ d^q 


^iqr 


2 / 8 3 2 , 2 

m^J OTT-^ 2 


By writing the exponential form of sin(...), we get 

.2 r ^3^ giq r 


27 4^2 j g ^3 q2 _|_ ^2 


T(q) 


By defining 


U(,K) := j 


d^q 


^—mR 


Svr^ + m? AttR 
with R = |R|, and by q —> —iV, we have 




4m^ 


T(-?V) 2 f/(r) - U{r+) - U{r_) 


with r-t = r ± A, and 

T(-iV) = W + (2S^ - 3)V^ - 2(S • V)^ + 6(p x S) • V 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


We mention that, for A = 0 the potential vanishes; this happens in the following cases: 
1) 0 = 0, 2) p = 0, and 3) p || 0. It is reasonable to see the behavior of potential for 
small noncommutativity parameter, defined here by A -C r and Am -Cl. In this limit, 
the first surviving terms are given by: 

^ 27 W = ^T(-*V)(A ■ VfU{T) + O(A^) (70) 

Recalling that for a function /(r), dif{r) = XiVrf, with = r~^dr, and using 

(p X S) ■ r = (r X p) ■ S = L ■ S, (71) 

V^U{r) =m‘^U{r)-6{r), (72) 
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with L as the total angular momentum, we get the expression for potential 


^27(r) = ^lm^{l + 2S^) 


V 


A2V,+ (A-r)V^V, 


+ 2(A . S)^ V.V. + (A ■ r)^(S ■ r) V.V.V.V. 


4(A ■ S) (A ■ r) (S ■ r) + A2(S ■ r)^ + ^^(A ■ r)' 
+6 A V, V, + (A ■ r) V, V,V J (L ■ S) 


\/ f \/ f' "p 


+2(p X S) ■ A (A ■ r) VrV,. 


t/(r) + D.D. + O(A^) 


(73) 


in which A = |A|, S = |S|, and D.D. is for the distributional derivatives of function U{r), 
containing 5-function and its derivatives; we calculate and present the explicit expression 
of D.D. in Appendix A. For sake of completeness, we just present the relevant expressions 


WrU 

VrVrU 

VrVrVrU 

VrVrVrVrU 


e mr + 1 
dvr ’ 

^-mr _|_ g 

dvr 

Q-mr _j_ g^2^2 _|_ _|_ ^5 

dvr r’^ 

e-mr + dSm^r^ + lOSmr + 105 

dvr 


(74) 


We make comments on the potential given by dZSl). First we mention that due to r’s in 
the inner products, the effective lowest power is r“®. Second, the strength of potential, 
through the definition of A, depends on momentum. Third, let us consider the spin- 
independent part of the potential, that is setting S' = 0: 




r = 


g2 g-mr 

d Att 


A 


2mr -|- 1 


+ (A 


-|- 3mr -|- 3 


(75) 


We mention that m = 0 limit of above expression is well defined. It is known that in 
noncommutative field theories particles behave as electric dipoles Uni CHI 1221123 123. 
The electric dipole depends on the strength of noncommutativity parameter as well as 
the momentum, and is perpendicular to both of them, and is given by d = ^eO x p. 
For the two-photon system, in center-of-mass frame, for which pi = —p 2 = p, we have 
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di = —d 2 = d. The potential for a system of two electric dipoles di and d 2 is given by 


hdipoles 


1 1 
dvr 


di ■ d2 


3(di ■ f)(d2 ■ f) 


(76) 


We see that two expressions m and dZED are eqnivalent for m = 0 and d = \eX. 
In fact the expression m is the potential of two anti-parallel dipoles in a theory in 
which the potential of a charged particle is given by the so-called Yukawa potential: 
^(^) = principle, one could justify that the potential is in fact that for 

two anti-parallel dipoles, included by spin-orbit and spin-dipole interactions in a Yukawa 
type theory. Finally, we mention: 


r = ^(p X 6 >)-r = ^6/-(r X p) = ^ 0-L (77) 

that can be inserted in the relevant parts of potential (HBD- It is the famous 6-L coupling, 
previously found in studies concerning the implications of noncommutativity in low energy 
phenomena UnillHlES. 


4 On Existence Of Bound States 

Having the effective potential, the starting point for studying the bound state problem is 
the Schrodinger-type equation by the Hamiltonian: 

H = 2m + H 2 h (78) 

in which m is the constituent mass, and H 2 h is sitting for the Hamiltonian capturing the 
dynamics of two-body system. For example, in the glue-ball case H 2 h usually consists 
three parts: the kinetic term, the potential term coming from perturbative calculation, 
like (HD), and string potential. The string potential usually is taken in the form V^tring = 
2 m(l -e-^’’), in which /3 is related to the tension of string stretched between the gluons. 
The formation of strings is expected from simulations on lattice, as well as the conhnement 
hypothesis EllllEI. Due to lack of analytical solutions, approximation methods, specially 
the variational method, appear to be practically useful 0111 ini. We mention that without 
any reliable estimation on the value of constituent mass, all efforts for the evaluation of 
bound state properties, such as mass and size, do not get any dehnitive result. There have 
been lots of theoretical and numerical efforts, like those done using the lattice version of 
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theory, together with phenomenological expectations, to estimate the mass of constituent 
gluons. 

Comparing to the case with glue-balls, the situation is more difficult in any study 
of photo-balls of noncommutative QED. First, by the present experimental data we just 
can suggest an upper limit for noncommutative effects, leaving unspecified 9. Second, at 
present neither we can say anything about the value of constituent mass, nor how it varies 
with other parameters, specially 6. In this sense, any study can not yield definitive result 
or suggestion for the quantities we like to know about photo-balls. 

Here we try to formulate the dynamics based on the effective potential obtained in 
previous section. Based on this formulation, we specially present a proof of existence for 
the bound states. Since the issue of possible formation of string-like objects in noncom¬ 
mutative QED is not in a conclusive situation, we do not consider a string potential in 
this work. We remind that by including the string potential the existence proof of bound 
states would be a trivial task. Also as the potential dZSD is very complicated for study 
the possible bound states, we restrict here ourselves to S' = 0 case; we also ignore D.D. 
terms. 

So we have potential and for sake of definiteness, we take the vector 6 in z 

direction, that is 0 = 6^ z. It is more convenient to work in cylindrical coordinates (p, cj), z), 
in which the kinetic energy, recalling that the effective mass in relative motion is m/2, is 
T = If (p 2 + p 202 

= + Py) = (P9) 

and also 

(A.y = l(A.y = f3mW.#'^ (80) 

in which r is the distance between two photons, r = ^+ z"^. We see, while the contri¬ 
bution coming from the velocity p always yields an attractive force, the contribution from 
angular velocity 0 depends on the ratio could be attractive or repulsive. In fact 

the ratio as represents how much the photons move off from the plane z = 0, also 

determines the relative orientation between r and the components of electric dipoles gen¬ 
erated due to velocity 0. We recall that the relative orientation of dipoles and the position 
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vector appears in dipole-dipole potential (|7U|). By these all we have the Lagrangian 


L = T -V = ^m[l + a/i(r)]p2 + [l + a(/i(r) - pV 2 (r))] 0 ^ (81) 


4 L ' /jr ■ 4 '4 

in which a = is a constant, and 

D47r ’ 

.mr + 1 


fiir) = e" 


, Ur) =-~ = e-”'’-- 
r or 


m?r’^ + 3mr + 3 


(82) 


We mention that the hrst two terms are positive dehnite, while third one can be negative, 
zero and positive. The coordinate 0 is cyclic, and hence its momentnm, given by 


or 1 
ocp Z 


(83) 


is a conserved qnantity, that we show with K; we see later that in qnantnm theory K 
should be an integer. One can hnd the effective theory for coordinates p and ;s, by 
eliminating 0 by using the Routhian R I 2 nj , as 


L pz R L 


1 1 
= -rn]! + afi{r)\p^ +-rnz^— 




4'"’^ ' ' d"'"’ mp2[l-h a(/i(r) - p 2 / 2 (r))] 

in which we recognize the potential 


(84) 


Ves{p,z) = 


mp^ [1 + a(/i(r) - p'^f 2 {r))] 
It is useful to mention the properties of Veg: 

• It goes to - 1 -cx) for p = 0 and 2 ; 7 ^ 0. 


(85) 


• It is 0 on p = 2 ; = 0. 

• It goes to ±cx) around the curve p(p, z) ■.= 1 + a[fi{r) — p‘^f 2 {r)) = 0. 

In Fig. 3 we have presented three plots of Wfr in pz-p\dja.e for m = a = 1, m = 10a = 10, 
and a = 10m = 10. We see that Ves goes to —00 and -|-cx) inside and outside regions 
defined by the curve g{p,z) = 0, respectively. We mention also, as the plots suggest, the 
dynamics on 2 ; = 0 plane is unstable; that is a small velocity i 7 ^ 0 hustles particles out 
of 2 ; = 0 plane. 
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Figure 3: Plots of Ves, for m = a = 1, m = 10a = 10, and a = 10m = 10. 


Before staring the discussion on quantum theory, let us have another look to the 
original Lagrangian (IHIl). We mention that the Lagrangian is in the form of a pure kinetic 
term, represented by means of a metric gij{x) as 

1 7Tl 

L =-—Qijix) x^x^, ( 86 ) 

in which x* = (p, (f), z), and 

9 ii{p, z) = 1 + a/i(r), P22(p, z) = [l + a(/i(r) - p^f2{r))] , 

933 = 1, 9ij = 0,^7^ j. (87) 

We remind that although the Lagrangian is looking like a pure kinetic term, since one of 
the components of metric, P 22 , changes sign, we can have negative energy states, among 
them there are bound states. By this interpretation of Lagrangian, the Hamiltonian of 
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quantum theory is simply gained: 


H2^.= 



1 

m 



VIdetsI g'^ 8j 


( 88 ) 


in which degg is the determinant of Qij. As Qij is diagonal, = 1/gij, for non-zero g^^^s. 
We mention deg g and components of g"^^ are independent of coordinate 0, and so we hnd 


^ 2 b = - 


1 

m 




+ 


a / 9n\g22\ 


A 


\/9ii\ 922\ dz 



(89) 


Using the separation of variables, we choose the wave-function 2:)<h(0), 

with <h(0) oc e*fo, and due to single valued-ness of wave-function, I should be an integer. 
So we replace —dl by —— in above, getting 


=- 


1 

m 


a/ fi'll|fl'22| 


-A 


k 


221 


9n 




+ 


a / gii\g22\ 


-rdz 


a/ gii\g22\ dz - \ (90) 


922 


in which means the Hamiltonian for states with specihed value for 1. By comparison, 
we see that the classical counterpart of integer number / is K. Also we mention that 
k/{mg 22 ), as expected, is sitting for Uefr in quantum theory. Now let us choose a trial- 
function f{p,z), that vanishes outside the curve g{p,z) = 0. We consider the quantity: 


ifWLlf) = 


— ^ 1 , 


inside g{p,z)=0 

fA2j =: E 


f*{p, z) (^hf 2 b/(P> ^)) VI detpl dp dz 


f 




(91) 


in which Ai j and A 2 J are two numbers independent of 1. We mention that, since f{p, z) 
vanishes in r — 00 , the contribution coming from the first two terms of iJ^b) faking into 
account the minus sign in front, is positive. The contribution from the last term of hf^b) 
reminding the dehnition of f{p,z), is negative. So for this kind of trial-function, Aij 
and A 2 J are positive. Here we make comment on the existence of bound states, at least 
for some ranges of /. We mention that for sufficiently large values of /, for a hxed trial- 
function /(p, z), Ef^i can be negative. In fact one can, by increasing /, lower Ef^i as much 
as wants. Now, by variational theorem we know that Ef^i is an upper limit for the lowest 
energy, and so we expect that for states with sufficient large /, there should be negative 
eigenvalues for Hamiltonian El\^. Showing these negative eigenvalues by E^^i, and the 
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corresponding eigenfunctions by 'ijjn,i{p, z), with n as for the possible quantum numbers, 
we have 





(92) 


with as the Laplacian in p; 2 -plane, given by 


V" = 


\/ 9ii\922\ 


--On 


w. 


22 


^11 


On 


+ 


\/ 9ll\922\ 


--dy 


\/ 9ii\ 922\ dz 


(93) 


Now, since outside of the curve g{p, z) = 0 the potential Ves is positive definite, the 
coefficient of in the right-hand-side is also positive. As r —> cx) belongs to the outside 
of the curve g{p,z) = 0, by the properties of spectrum of V^, we expect 'ipn,i\r^oo 0, that 
is is representing a bound state. Physically we expect that for the negative eigenvalues, 
the wave-function should be localized along the well inside the curve g{p,z) = 0, as 5^22 
approaching zero from below. 

The manner we approved the existence of bound states can be used, by increasing I, 
for reasoning that there is no lowest energy state: the eigenvalues are unbounded from 
below. We recall that the potential dZSl) is obtained under the assumption that A -C r. As 
A = I p X 0, we see that for large values of momentum, A may be comparable, and even 
bigger than r. One situation that might invalidate the assumption A -C r can happen for 
very large values of I, corresponding to large value of K in classical theory. In such cases 
one should consider the original potential (EH). We remind that, although the absolute 
least energy is meaningless to be found under the approximation A -C r, the least value 
of energy is still meaningful for states with specihed value for 1. 

The other issue is about states with eigenvalues bigger than the maximum of potential 
inside the curve g{p, z) = 0. We mention that an infinite tall wall has surrounded the inside 
region, and the question is if the wall can make the possibility for forming bound states. 
In fact since the thickness of wall behaves like 1/h, with h as height, by considerations 
coming from the WKB approximation for tunnelling effect, we expect that the particles 
with positive energies can escape from the inside region. This situation is similar to the 
situation in one dimensional problem with potential V{x) = l/{x — Xq), for which by 
WKB method one finds a finite probability expression for tunnelling of positive energy 
particles. 

As the hnal point, we make comment on the possible values of spin and 1. The state of 
a two-photon system should be symmetric under the exchange of photons. A two-photon 
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system can have 0, 1 and 2 as total spins, as for the hrst and last ones the spin states 
are symmetric, and for the second is anti-symmetric. Here the exchange of two photons 
means z ^ —z and 0 —> 0 -f tt. By considering the spatial dependence of wave-fnnction, 
we have the followings for allowed spins and l\ 

5 = 0,2, / = 0,2,4,-■■ 

5=1, / = 1,3,5,--- (94) 

5 Conclusion And Discussion 

We mention that the transformations of gange held as well as the held strength in a non- 
commntative space look like to those of non-Abelian gange theories. Besides we see that 
the action of noncommntative QED contains terms which are responsible for interaction 
between photons, again as the sitnation we have in non-Abelian gange theories. There 
is another observation that promotes the formal similarities of noncommntative and non- 
Abelian theories to their behaviors, that is the negative sign of 0-fnnction, which manifests 
that these theories are asymptotically free HSlIIll. The above mentioned observations 
make it reasonable to stndy whether and how the photons of noncommntative QED can 
make bonnd states. Also these observations make it reasonable to see if the techniqnes 
developed for QCD pnrposes can also be nsed for noncommntative QED. Here we nsed the 
so-called potential model, developed on the constitnent glnon pictnre of QCD glne-balls. 
The basic ingredient of potential model is that the self-interacting massless gange particles 
may get mass by inclnsion non-pertnrbative effects. By calcnlating the amplitnde for the 
scattering process between two massive photons, we extract the effective potential that 
is expected captnre the dynamics of constitnent photons. Using this effective potential, 
we formulate the Hamiltonian dynamics, by which arguments are presented in favor of 
existence of photon bound states. 

As possible photo-balls, like their glue-ball cousins, are non-pertnrbative in nature, 
it is expected that lattice version of noncommntative QED should appear as one of the 
natural ways to study photo-ball’s properties. It is remarkable to remind that ordinary 
QED on lattice develops an area law, suggesting a stringy picture for force, for two charged 
particles EJI- There are suggestions for lattice version of noncommntative gauge theories 
PH] • Specially, the finite N version of the theory is promising for numerical and simulation 
purposes. Recently, there have been a few works reporting the preliminaries results by the 
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lattice version of theories [2^1 • There are other suggestions for non-perturbative dehnition 
of noncommutative QED |d()j . 

By the current experiments there has not been any signal for possible noncommutativ¬ 
ity. So the common expectation is that the evidence for noncommutativity, if any, should 
modifies the processes that occur in energies much higher than those presently available. 
It is why that by present experimental data one can just suggest an upper limit for 
noncommutative effects. There has been another suggestion that the noncommutativity 
effects may appear due to applying sufficiently strong magnetic held on samples contain¬ 
ing moving charged particles. It would be extremely interesting if noncommutative view 
let us know something new about relevant phenomena |dlj . 

Acknowledgement: A. H. F. is grateful to M. Khorrami for very helpful discussions 
on the distributional derivatives, and also for extremely useful discussions on the bound 
state problem. 

A Distributional Derivatives 

Here we calculate the distributional derivatives. One very helpful reference is [21] • First 
we consider 

p-mr 

d,d, - (95) 

r 

The distributional derivative can be calculated by its effect on a test-function 0(r) 

/ p-iTT-r , /p-mr , r p-mr 

= lim j - -djdj'0(r)dV (96) 

e—>0+ JT 

in which d^r = r^drdO. The limit above does exist because the integral in the second 
line, due to in dV, is hnite. r = 0 is excluded from the last integral, and so we can use 
the identity 

p— mr >p— mr ^ p— mr 

——didj(j) = di(^——dj(j)J-di——dj(j) (97) 
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by which we have for lij = (^didj^-^, 0 

.^-mr . r -mr 

di ^-(9j0 j dV — j di - dj(j) dV 

nr . p ^—mr 

— dj4>jd\— J di - dj 0 d V 


Ia = lim 

£^0+ 


L J r>e 


= lim 

e^0+ 

= lim 

e^0+ 

= lim 

e^0+ 


V- e 


L J r>£ 


\_j r=e ^ 


-dj(j){—ei ■ f)r^df2 


' r>£ 


di -(9j'0dV 


r 


p—mr n „—mr 

-0,'0 njr^dhl — / di -0^0 dV 

^ J r>£ 


r 


(98) 


in which dj is for unit vector in x, y and directions, and f = (ni,n 2 ,n 3 ). The hrst 
integral in last line is proportional to e and so vanishes in the limit. So we get 


In = — lim 


‘■IJ 

£—>0+ J r>£ 

By repeating the steps above, we arrive at: 


di - djcj) dV 


(99) 


In = lim 

£^0+ 


= lim 

£^0+ 


(t)di- 


riir^dVL + 


(t)didj- 


dV 


' r>€ 


— (pie) J + J (pdidj -dV 


( 100 ) 


for getting it we used dif{r) = riidrf^r), and also keeping the integrand of hrst integral 
to the hrst non-vanishing order in e. We know 


niUjdfl = A 25 ij 


( 101 ) 


and A 2 can be calculated by taking trace of the both sides, yielding A 2 = By these 
all we have 


dvr /* G 

lij = — 5-0(O)5jj + lim 0 didj — 




dV 


( 102 ) 


r>€ 


The limit above exists, by using the fact that the value of function in origin is constant 
and independent from the hi = Q{6, cp), and recalling 


{Sriifij — (5jj)dr2 = 0 


(103) 
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One can remove the limit by respecting the order of integrations. By these all we get 


didj- 


-y(5ij(5(r) +pf 


didj- 


( 104 ) 


in which “ pf ” is sitting for pseudo-function, that here simply means that in integrals the 
integration on solid-angle should be done before radial one. 

Now we come to 


Hjk 


(didjdk 





(9i(9j(9fc0(r)dV 


I P~ 

— lim / -(9i(9j(9fc0(r)dV 


(105) 


in which again the limit exists because the integral in the second line is hnite. By repeating 
the steps for calculation lij, and the replacement 


0 ( 5 ) = 0(0) -|- e f ■ V0(O) -|- O(e^) 


(106) 


and using 


/ 


Tl{TljTlf^Tlp(l.^ ^4 ^ij^kp H” ^ik^jp ^ip^jk 


/ 


njdO 




( 107 ) 


with A 4 = ^, one reaches to 


Hjk 


dTT 

T 


Sijdk(l>{0) -F 6jkdi(l){0) + (5fci(9j0(O) 


I 

-I- lim / (j) didjdk -dV 

^^o+Jr>e r 


(108) 


for which again the limit above exists while one is careful that the integration on dfl 
should be done hrstly. By these all we get 


didjdk- 


An 

'T 


5ijdk + Sjkdi + Skidj (5(r) -F pf 


djdjdk- 


(109) 


in which again “ pf ” simply means in integrals the integration on solid-angle should be 
done before radial one. 
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Now we come to 


hjki = {didjdkdi—^,(t)y.= (^—^,didjdk 


( 9 i( 9 j( 9 fc( 9 i 0 (r)dV 


= lim 


didjdkdi(t>{T:)i\ 


( 110 ) 


r>e 


in which again the limit exists because the integral in the second line is hnite. By repeating 
the steps for calculation lij and lijk, and the replacement 

^2 


and using 


(P{e) = 0(0) + er- V0(O) + -{r ■ V)20(O) + 0{e^) 


J riiUjnknpnqdVL = 0 

TliTljTl}^TliTlpTlq^^ -A.Q^SijS}^lSpq ~\- ^ik^jl^pq ~\~ ^il^kj^pq 

~\~5ij5kp^lq H“ ^ij^kq^lp H“ ^ik^jp^lq 
jq^lp H“ ^il^kp^jq H“ ^il^kq^jp 

+ 5ip5 kl^jq H“ SipSkjSlq “h Sip SkqSjl 

+ 6iq6 klSjp Siq ^kj^lp ^iq^kp^jl^ 


( 111 ) 


( 112 ) 


with Aq = 


one reaches to 
Att 


lijkl - - 


15 

207r 


^ij^kl A~ Sik^jl “ 1 “ ^il^jk 0 ( 0 ) 


iOtt r 


51 


^ij^kl A~ Sik^jl SiiSjk V 0(0) 


Sijdkdi + Sikdjdi + Siidjdk + hAdj + 5 jididk + ^kjdidi 0(0) 


lim 

e^0+ 


0 didjdkdi -dV 


(113) 


r>£ 


for which again the limit above exists while one is careful that the integration on dfl 
should be done firstly. By these all we get 


didjdkdi 


Ati 

’l5 

207r 


1 Ott 

^ij^kl “ 1 “ ^ik^jl A~ Sildjk <^(r) p, .. I ^ij^kl ^ik^jl “ 1 “ ^il^jk I V *^(^) 


51 


51 L 


Sijdkdi + Sikdjdi + Siidjdk + SkAdj + SjidA + hjdidi ( 5 (r) 


+pf 


didjdkdi- 


(114) 


28 



in which again “ pf ” simply means in integrals the integration on solid-angle should be 
done before radial one. The combination V'^didj is simply S^^'didjdkdi. 
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